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^ ■ Abstract 

V^ . We present an analytic proof of the integrals for astrophysical thermonu- 

O ! clear functions which are derived on the basis of Boltzmann-Gibbs statistical 

Q I mechanics. Among the four different cases of astrophysical thermonuclear 

"^ I functions, those with a depleted high-energy tail and a cut-off at high ener- 

CL(' gies find a natural interpretation in q-statistics. 

-(— > ■ 

g ■ 1 Introduction 

.5^ I One of the first applications of Gamow's theory of quantum mechanical po- 

/\ • tential barrier penetration to other than his analysis of alpha radioactivity 

j^ ■ was in the field of thermonuclear astrophysics [1,2]. Atkinson and Houter- 

mans proposed that the source of energy released by stars lay in thermonu- 
clear reactions taking place near their centres where the motion of nuclei was 
supposed to be in thermal equilibrium. The state of hot stellar plasmas is 
such that only the lightest chemical elements could contribute because of the 
Coulomb repulsion between nuclei. In effect, the rate of energy production is 
goverend by the average of the Gamow penetration factor over the Maxwell- 
Boltzmann velocity distribution [3]. The thermonuclear reaction rate is the 
coefficient in the rate equation that is used to describe the change of chem- 
ical composition of hot plasmas. Similarly fundamental, the inverse of this 
coefficient is the characteristic time scale for the respective thermonuclear 



reaction. Nuclear and neutrino astrophysics are vibrant fields with many 
experiments collecting data continuing exploring energy production in stars, 
particularly the Sun, deduced from neutrino measurments [4]. 

Commonly, hot stellar fusion plasmas are described in terms of Boltzmann- 
Gibbs statistical mechanics based on the entropy Sbg = ^kYl^iPdnpi. 
Thermo dynamical and ststistical description of nonextensive systems may 
require a generalization of Boltzmann-Gibbs thermostatistics. Examples of 
physical systems or processes where Boltzmann-Gibbs approach seems to be 
inadequate are large self-gravitating systems and hot and turbulent plasmas. 

An ultimate generalization of Boltzmann-Gibbs thermostatistics is due 
to Tsallis [5] based on the following expression for entropy 



9 - i=i 



S„ = ^^^^f^, qeR,Si = SBG, Ep^ = 1' (1) 



with the following expression for equal probabilities 

14/1-9 - 1 
Sg{p, = 1/W, V,) = k—^ = klugW, (2) 

S,{p^ = l/W,y,) = klnW = SBG. 
The optimization of the entropic form (1) under the restriction 

< H >.- ^0^ = U„ (3) 

2^j=i Pj 

where < . . . >i=< . . . > pi / Y^=iP'j is the escort distribution, {Ei} are the 
eigenvalues of the Hamiltonian H with appropriate boundary conditions, and 
the internal energy L'^ is a finite fixed value, yields 

-(5^{Ei-Uq) 

P. = ^^ , (4) 
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with 

w 



^.^Ee/'^^^-'''^ /3. 
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"<? ■, Hq — ^w <1 ' 

i=i ^i=i t^j 



(3 = 1/kT. Tsallis [5] verifies that q = 1 recovers Boltzmann-Gibbs weight, 
g > 1 implies a power-law tail at high values of Ei, and g < 1 implies a 



cut-off at high values of Ei. Subsequently, this behavior is reflected in the 
equilibrium distribution of energies. In the following, a method for obtain- 
ing closed-form representations of thermonuclear reaction rates (astrophys- 
ical thermonuclear functions) is developed, taking into account the above 
described behavior of q-statistics, these are the cases g > 1 (dubbed "de- 
pleted") and g < 1 (dubbed "cut-off") for the energy distribution. Thus, 
q-statistics provides a natural physical interpretation for hot plasmas where 
deviations from the Maxwell-Boltzmann distribution are expected. 

2 Definition of astrophysical thermonuclear 
functions /^(z/ — 1, a, ^, p) 

By employing a statistical technique, Haubold and Mathai [6] have estab- 
lished the following integral formula for the derivation of closed-form repre- 
sentations for four astrophysical thermonuclear functions: 

/ y^-^ exp[-ay-zy-f]dy=^I'^{iy-l,a,z,p) (5) 
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{u+r+1,1) 



^ l(0,l/p)(i/+r,l) 



(Pt^ 0),/or <i< oo; 



= -^0,2 [az^^''\{o,i/p),{u,i) (p 7^ 0), for d = oo, 

where Reiy) > 0, Re(a) > 0, Re(z) > 0, Re(p) > 0. It has been shown 
that by the application of the Mellin-Barnes integral rerpresentation of the 
exponential function, the given integral can be transformed into a Mellin- 
Barnes type integral representing an H-function [7]. We also deduce the 
value of this integral when p < 0. Results for all the four astrophysical 
thermonuclear functions are also obtained. The integral I^{z,t,v) [8, see 
also 2,9] is evaluated in a generalized form. For the sake of simplicity and 
continuation we adopt the same notations as used in [6, see also 2,9]. In order 
to prove the integral (5), we employ the following Mellin-Barnes integral 
representation for the exponential function, namely [10,11] 

exp[—x\ = - — / T{s)x~^ds, \x\ < cxd. (6) 

ZTIUJ Jl 



where L is a suitable contour {u = (—1)2), interchange the order of inte- 
gration, the given integral, denoted by /'^(z/ — 1, a, z, p)), transforms into the 
form 

I^u - 1, a, z, p) = -^ / Tis)z-^ /' y''+^-'-'expi-ay)dyds (7) 

Zttu Jl Jo 



We know that 
rd 
/o 



I'd 

/ exp{-ay) y^'^dy 
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(8) 



^^T{u + r + l){r)\ ^ ^ 

^(^/)/a^ for ci = 00 i?e(a) > 0, Re{u) > 0. 

By virtue of the above values of the y-integral for d < 00 and for c? = 00, it 
is found that 



I'^{u-l,a,z,p) 



d^ ^ i-ady r r{s/p)r{u + s + r)z-^/P ^^ 

(p 7^ 0) ioi d < 00 

^^ / r(-)r(z/ + s)a-^2-^/''rfs, 
2T[ujp Jl p 

(p 7^ 0) for (i = 00 



The integral formula (5) now readily follows from the equations (9), if 
we interpret the above contour integrals in terms of the H-function, which is 
defined by means of a Mellin-Barnes type integral in the following manner 

[7,p.2]. 
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where 



e(0 



n"L,r(6, + 5,0 n"=,r(i - a, - a,o 



TT'' 

jr'=m+l 



^T{l-h,-B,i) m^^,V{a, + A^i) 
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an empty product is interpreted as unity ; 171,11, p,q G A^o^ith < n < 

p,l < in < q, 

Ai, Bj G R+, ai, bj G -R or C{i = 1, . . . ,p;j = 1, . . . ,q), such that 

Ai{bj + k) ^ Bj{ai-l- l){k,l e N - 0;i = 1, . . . ,n;j = 1, . . . ,m). (12) 

We employ the usual notations: A^o = (0,1,2...); R = (—00,00),/?+ = 
(0, 00) and C being the complex number field. L is one of the contours de- 
scribed in the monograph by Mathai [11]. In case we make the transformation 
p = —1] with 1] > in {9), then the following result holds. 



I'^{u-l,a,z,-ri) 
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r=0 



A/ri 



(l-u-r,l) 



(13) 



1] ^ 0, for z < 00; 
"-"1,1 



A/v 



1] 
1] j^ 0, for 2; = 00, 






(14) 



where Re{u) > 0,Re{a) > 0,1] > 0,Re{z) > 0. When p is real and rational 
then the H-functions appearing in (5), (13) and (14) can be reduced to G- 
functions by the application of the well-known multiplication formula for 
gamma functions [12]: 



r(inz) 
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1,2,... (15) 



3 The astrophysical thermonuclear functions 

Ii{z,iy) and l2{z,d,jy) 



As an application of the results (5), (13), and (14) we deduce the values of the 
following two thermonuclear functions Ji [8] and I2 ["cut-off", 8,13,14,17], in 
terms of which the astrophysical thermonuclear functions are expressed [8, 
see also 2,9]: 



h{z,u) =^ J y" ^exp[-y Y/^]dy 



(16) 
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where Re{z) > 0,Re{u) > and GoslOis the Meijer's G-function [10]. 



h{z,d,v) 



def 



71 



y" ^exp[~y - zy ^''^]dy 

T/2 2^ ^rvr'-^i,3 



1^0 (0! 



Ad 



[v+r+l) 
(i^+r), 0,1/2 



(17) 



where d > 0, Re{z) > 0, Re{u) > 1. 



4 The astrophysical thermonuclear functions 

Is{z, t, jy, /J,) and I^^z, S, 6, u) 

We now show that [81 



h{z,t,u,^) 



def 



y" 'exp[-{y + z{y + t) ^'}]dy 

oo j.r 



(18) 



il+u,fi),{u+r+l,fj.) 
(0,1),(1+1',/^),(1,m) 



r=0 



r ! 



(l-'':/i).{'^.M) 



where Re{v) > 0, Re{z) > Oand /i > 0. To prove (18) we observe that in 
view of the formula (6), the value of the integral is equal to 



1 r r°^ 

—- / T{s)z-' / y^-^[exp{-y)]{y + trdyds 
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(19) 



Evaluating the y-integral in terms of the Whittaker function, defined in [12, 
p.255] and [12, p.257]: 



\E'(a, c; z) 



1 rOD 

-- / x^-^e-'^'il + xy-^-^dx, Re{a) > 0, Re{z) > 0. (20) 
(a) JO 

^^^ " ""^ <l>fa, c, 2) + ^M^z^-'<!>{a - c+ 1,2 -c;z) 



T(a) 



r(a-c+l 

r(l-c) ^ (a),z'' r(c-l) ,_^^ (g-c+l),^'- 
r(a-c+l),^o(c)r(r)!+ r(a) "^ .^^ (2 - c).(r)! ' 



where c is not an integer and <l>(a, c; z) is Kummer's confluent hypergeometric 
function [12, p. 248], the integral (19) reduces to an elegant formula 



2noj 



I r(s)^(z/, z/ + s/i + 1; t)z-Vds. 



(21) 



The integral (21) can be evaluated by expressing the Whittaker function 
appearing in its integrand in terms of its equivalent series given by (20) 
above, and reversing the order of integration and summation. Thus the 
given integral finally transforms into the form: 



h [z, t,iy,fi)=t 2^ 

r=0 

oo 



T{s)T{-u - s/i)r(z/ + sfi + ^) ^-s^sf,^^ 



+ E 



f 



r=0 



ir)\ 271U Jl 



r)! 2-KUj Jl r(— s/x)r(z/ + r + s/i + 1 

r(s)r(z/ + Syu)r(r — syu)r(i — v — sfi) 



r(— s/i)r(i + r — u — sfi) 



-z '^ds. 



(22) 



which on interpreting with the help of the definition of the H-function (10) 
establishes the desired result (18). If we set /i = 1/2 in (18), then by an 
appeal to the duplication formula for the gamma function (this is equation 
(15) for m = 2), the H-functions reduce to G-functions [10] and consequently, 
we obtain 



h{z,t,u) 



def 
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(23) 



where Re{v) > 0, Re{z) > 0. Next we will prove the formula ["depleted", 8, 15, 16, 17] 



Ii{z,6,b,u) 



def 



E' 

r=0 



y'^ ^exp[— ly + by^ + zy2\]dy 
i-b/z) 



-H, 



r ! 



2,0 
0,2 



(r,l),(j/+r5+r/2,l/2) 



(24) 



(25) 



where Rei^u) > 0, Re{z) > 0, Re{b) > and 5 > 0. To estabhsh the integral 
formula (25), we see that in view of (6), it can be written as 

/•oo 1 f 

h{z, 6, b,u)= y'^'^expi-y)—- / T{s)y^/\by'+'/' + zY'dsdy (26) 
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On employing the formula 



~ ia) 



(l + ^)"" = E73[(-^r, kl<l, (27) 

r=o \ /• 

and reversing the order of integration and summation, the equation (26) 
transforms into the form 



oo 



which, when interpreted with the help of (10), yields the desired result (25). 
It may be noted that the result (25) can be expressed in terms of the G- 
function in the form : 



\r 



1 - (-26/^ 3,0 
7rV2 ^^ (r)\ °'^ 



hiz, s,b,iy) = —-^J2 ^7:;x]— ^0 



r=0 



Z^ 



r r + 1 



(iy+r5+§) 



(28) 



where Reiy) > 0, Re(b) > 0, Re{z) > and S > 0. Further it is interesting 
to observe that, as b tends to zero, (28) reduces to (16). 

5 Relation of astrophysical thermonuclear func- 
tions to Kraetzel functions 

It is not out of place to mention that the function represented by the integral 
(5) in case ;i; = oo has been studied by Kraetzel and is known as Kratzel 
function in the literature. It may be noted that Kraetzel [18] introduced the 
integral transform 

{K^J){x) = j^ Z;{xt)f{t)dt {x>0), (29) 

involving the kernel function 

Zl{x) = / f-^expl-ff - -]dt {p>0,uE C). (30) 

JO t 



Further he has also investigated the asymptotic behavior of Zp{x) as a; — > 
and a; — > cxD in [18]. 

Finally it may be mentioned here that the integral described here for 
z = oo plays not only an important role in the study of astrophysical ther- 
monuclear functions but also for Bessel-type fractional ordinary and partial 
differential equations [19]. 
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